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ABSTRACT 

We report the discovery of a new detached, double white dwarf (WD) system 
with an orbital period of 39.8 min. We targeted SDSS J163030. 58+423305.8 (here- 
after J1630) as part of our radial velocity program to search for companions around 
low-mass WDs using the 6.5m MMT. We detect peak-to-peak radial velocity varia- 
tions of 576 km s _1 . The mass function and optical photometry rule out main-sequence 
companions. In addition, no milli-second pulsar companions are detected in radio ob- 
servations. Thus the invisible companion is most likely another white dwarf. Unlike 
the other 39 min binary SDSS J010657. 39— 100003.3, follow-up high speed photometric 
observations of J1630 obtained at the McDonald 2.1m telescope do not show signifi- 
cant ellipsoidal variations, indicating a higher primary mass and smaller radius. The 
absence of eclipses constrain the inclination angle to i ^ 82°. J1630 contains a pair 
of WDs, 0.3Mq primary + ^0.3Mq invisible secondary, at a separation of ^0.32Rq. 
The two WDs will merge in less than 31 Myr. Depending on the core composition of 
the companion, the merger will form either a single core-He burning subdwarf star or 
a rapidly rotating massive WD. The gravitational wave strain from J1630 is detectable 
by instruments like the Laser Interferometer Space Antenna (LISA) within the first 
year of operation. 

Key words: binaries: close - - white dwarfs stars: individual (SDSS 
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1 INTRODUCTION 

Radial velocity surveys of WDs are the best way to iden- 
tify short period systems that may merge within a Hubble 
time and produce normal Type la or underluminous ".la" 
supernovae ( SNe. IWebbinkl 1 1984 llben fc Tutuko^l 1 19841 ; 



iBildsten et al.|[2007 'l. However, even large spectroscopic sur- 
veys like the SN la Progenitor Survey (SPY) have lim- 
ited s u ccess in identifyi n g mer ger systems (|Napiwotzki et all 
l200ll : iNelemans et alj l2005l h Many of the radial ve- 
locity variables in SP Y and othe r surveys turn out 
to be low-mass WDs l|Marshl Il995l : iMaxted et ail I2OO0J : 
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iNapiwotzki et alj|2007t ). Only 24 double WD systems were 
known in 2 005, and all but three have primary WDs with 
M ^O.5M l|Nelemans et aT]|2005h . 

The Galaxy is not old enough to produce low-mass WDs 
through sin gle star evolution, but they can form in unu sual 
systems (see lKilic et~al]|2007bl : INelemans fc Taurislll998l . for 
a review). The majority of low-ma ss WDs are found in bi- 
nary systems (I Brown et al.l l2011aT ) , where the progenitor 
main-sequence stars experience enhanced mass-loss during 
one or two common-envelope phases. Hence, low-mass WDs 
are prime targets for discovering short period merger sys- 
tem^ 

iKilic et al.l (|20ld . l2011ah and lBrown et all (|2010D have 
established a radial velocity program, the ELM Survey, to 
search for companions around known extremely low-mass 
(M ~0.2M Q ) WDs in the SDSS Data Rel ease 7 foot- 
print and the MMT Hypervelocity Star Survey (|Brown et al.l 
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l2006h . So far, the ELM Survey has identified 12 WD merger 
system s, tripling the number of systems known (jKilic et al.l 
l2011al ). Depending on the unknown mass ratios, these sys- 
tems may be the progenitors of stable mass-transfer AM 
CVn binaries, single helium-rich subdwarfs, extreme helium 
stars, or underluminous SNe. The ELM Survey is also re- 
sponsible for findi ng the first tidally distorted WD J0106 
jKilic et al.ll2011bh and the 12 min or bital period detached , 
eclipsing binary WD system J0651 (jBrown et al.l l2011ch . 
These two systems, J0106 and J0651, are currently the only 
known detached binary WD systems with periods shorter 
than 1 hr. All other known binary WD systems with P < 1 
hr are in interacting AM CVn systems. 

Here we present the discovery of a new detached bi- 
nary WD system with P < 1 hr found in the ELM Sur- 
vey. J1630 was re cognized as a low-mas s WD in the SDSS 
DR4 WD catalog i|Eisenstein et al.ll2006h . Our follow-up ra- 
dial velocity observations demonstrate that J1630 contains 
a pair of WDs with an orbital period of only 39.8 minutes. 
In Section 2 we describe our spectroscopic and photometric 
observations. In Sections 3 and 4 we constrain the physical 
parameters of the binary and discuss the nature and future 
evolution of J1630. We summarize our conclusions in Section 
5. 



2 OBSERVATIONS 

We used the 6.5m MMT with the Blue Channel spectro- 
graph to obtain medium resolution spectroscopy of J1630 
on UT 2011 March 1 and April 22-23. We operate the spec- 
trograph with the 832 line mm -1 grating in second order, 
providing wavelength coverage from 3600 A to 4500 A and 
a spectral resolution of 1.2 A. J1630 has g = 19.0 mag. We 
started our observations with 400 s exposures, but reduced 
the exposure time to 300 s after detecting significant velocity 
variations in back-to-back exposures. We obtain all observa- 
tions at the parallactic angle, with a comparison lamp ex- 
posure paired with every observation. We flux-calibrate us - 
ing blue spectrophotometric standards (|Massev et al.ll 19811) . 
and we measure radial velocities using the entire spectrum 
in the range 3700-4430 A and the cross-correlation package 
RVSAO. T he details of our d ata reduction procedures are 
discussed in lKilic et all (^OIOT ). 

Realizing that J1630 has a period similar to the tidally 
distorted WD J0106, we acquired high speed photomet- 
ric observations of J1630 to search for eclipses and el- 
lipsoidal variations. We used the McDonald 2.1m Otto 
Struve Telescope with th e Argos frame transfer camera 
jMukadam fc Natherl [20051 1 over 1 night in 2011 May and 

3 nights in 2011 July. We obtained 1819 images of J1630 
with a 1mm BG40 filter and 15 or 30 s exposures for a total 
exposure time of 9.6 hr. 



3 RESULTS 

3.1 The Orbital Period 

Table 1 lists our radial velocity measurements for J1630. 
We compute best - fit orb ital elements using the code of 
iKenvon fc Garcial (|l986T ). To verify the uncertainty esti- 
mates, we perform a Monte Carlo analysis using 10000 sets 



Table 1. Radial Velocity Measurements for J 1630 
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Figure 1. The heliocentric radial velocities of J1630 observed 
over three nights in 2011 March and April (left panels). The right 
panel shows all of these data points phased with the best-fit pe- 
riod. The solid line represents the best-fit model for a circular 
orbit with a period of 39.8 min and K = 288.1 km s . 

of simulated radial velocities (see iBrown et al.ll2010l ). We 
adopt the inter-quartile range in the period and orbital ele- 
ments as the uncertainty. 

J 1630 exhibits radial velocity variations with a semi- 
amplitude of K = 288.1 ± 4.9 km s~ and orbital period 
of P = 0.027659 ± 0.000043 d, or 39.8 ±0.1 min. Figure 
1 shows the best-fit orbit compared to the observed radial 
velocities. Due to our relatively long exposure times, our 
data do not sample the velocity curve well near maximum 
and minimum. Thus, we underestimate the velocity semi- 
amplitude. To correct for this underestimate, we sample a 
sine curve at the exact 21 phases of our observations with 



P/8 long integrations. We recover the exact period, but K is 
systematically underestimated by 2.7%. Thus, the corrected 
velocity semi-amplitude is K = 295.9 km s _1 . With this cor- 
rection, J1630 has a mass function of / = 0.07423 ± 0.00369 
Mq. The systemic velocity (after subtracting the gravita- 
tional redshift of 7.6 km s _1 ) is —18.7 ± 3.8 km s _1 and the 
time of spectroscopic conjunction is HJD 2455622.024130 ± 
0.000097 d. 



3.2 The Physical Parameters of the Binary 

lEisenstein et al.1 l |2006l ) derive T eff = 14850 ± 360 K and 
log<? = 6.89 ± 0.13 from t he low-quality SDSS spectrum of 
J1630. iKilic et al.1 <|2007ah find T cff = 14180 ± 940 K and 
loggr = 7.08 ± 0.07 based on a re-analysis of the same spec- 
trum. Our higher resolution and higher signal-to-noise MMT 
spectra provide a better measurement of the atmospheric pa- 
rameters of the visible WD. We perform stellar atmosphere 
model fits using synthetic DA WD spectra kindly provided 
by D. Koester. 

Figure 2 shows the composite spectrum and our fits 
using the entire spectrum (top panel) and also using only 
the Balmer line profiles (bottom left panel). To derive ro- 
bust statistical errors, we also perform fits to the individual 
spectra. We find a best-fit solution of T cff = 14670 ± 320 K 
and log g = 7.05±0.08 from the composite spectrum. Fitting 
only the Balmer lines, we obtain T e g = 15560 ± 590 K and 
loggr = 6.93 ± 0.05. These results are consistent with each 
other and the previous T e ff and logg estimates within the 
errors. The best-fit model using the continuum shape pro- 
vides a slightly better fit to the SDSS photometry (Figure 
2) and we adopt that model for the remainder of the paper. 
We n ote that all of the stu dies described above ar e done 
using iKoested (|2010h models. iGianninas et al l l|2011i ) find a 
systematic offset of log g « 0.1 dex betwe en their analysis 
and the analysis done using iKoesteri (|2010h models. Hence, 
the systematic errors in our logg measurement are ~ 0.1 
dex. 

Based on th e improved IPanei et al.l ()2007l ) tracks (see 
IKilic et al.ll2010h for ELM WDs, the visible WD in J1630 is a 
160 Myr old 0.30 ± 0.02 Mq WD with a radius R = 0.025 ± 
0.003 Rq. Its absolute magnitude M g = 9.8 corresponds to 
a distance of 0. 7 kpc. Based on five epochs from the USNO-B 
and the SDSS,[m unn et al.1 <|2004D measure a proper motion 
of (fi a cosS, /if) = (2.3 ± 3.5, —7.3 ± 3.5) mas yr _1 . J1630 is 
0.5 kpc below the Galactic plane and it has U = 29±11, V = 
— 13 ± 9, and W = — 11 ± 9 km s -1 with respect to the local 
standard of rest l|Hogg et al.ll2005h . J1630 is a disk star. 

The mass function implies a companion mass of 0.30 
Mq . Such a main-sequence companion would have been de- 
tected in the SDSS photometry. Therefore, the companion 
is a compact object. Based on the mass function alone, the 
probability of a neutron star (1.4-3.0 Mq) companion is 5%. 
However, n o milli-second pulsar companions are detected in 
the radio (Ag iieros et al.l [20091 ). Hence, the companion is 
very likely another WD. 

3.3 The Light Curve 

Figure 3 shows the fourier transform (FT) of the J1630 light 
curve. The light curve is affected by low frequency sky noise 
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Figure 2. Model fits (red lines) to the composite spectrum 
(jagged line, top panel) and to the Balmer line profiles of J1630 
(bottom left panel). The spectral energy distribution of J1630 
(filled circles, bottom right panel) compared to the best-fit model 
using the continuum shape (solid line) and using only the Balmer 
line profiles (dotted line). 

due to the changes in the atmosphere and the color differ- 
ences between the WD and the reference stars. The middle 
panel shows the FT pre-whitened by the main peak in the 
top panel. There are many aliases, but the highest peak is 
at the orbital period (dashed line). The FT for the brightest 
comparison star (bottom panel) does not show any signifi- 
cant variations at the same frequency. 

Figure 4 shows the Argos light curve of J1630 (top 
panel) folded over the best-fit orbital period of 39.8 min. The 
bottom panel shows the same light curve binned into 100 
orbital phase bins. Unlike J0106 and J0651, J1630 does not 
show any significant ellipsoidal variations. The amplitude 
of the ellipsoidal effect is roughly 5f e u = (m,2/mi)(ri/a) 3 , 
where a is th e orbital semi-majo r axis and ri is the ra dius of 
the primary (jZucker et al.ll2007l : Isirporer et al.ll2010h . Even 
though J0106 and J1630 have almost the same orbital pe- 
riod, J 1630 is about a factor of two more massive and a fac- 
tor of two smaller in size compared to J0106. Therefore, the 
expected amplitude of the ellipsoidal variations is <0.1%, 
consistent with the observations. 

We expect to see «0.24% amplitude sinusoidal varia- 
tions due to the relativistic beaming effect. Fitting a sinu- 
soid, the lightcurve is best explained by 0.17% ± 0.05% vari- 
ations at the orbital period (see the bottom panel in Figure 
3). These variations are consistent with the expected beam- 
ing signal. To verify that the observed variations in the light 
curve are not a statistical fluctuation, we perform a boot- 
strap analysis. We create 10000 synthetic light curves by 
random sampling with replacement from the observed light 
curve. This analysis shows that the probability of measur- 
ing a random signal with ^0.17% amplitude at the orbital 
period is 2%. Based on Figures 3 and 4 and our statisti- 



M. Kilic et al. 



0.4 



0.4 




0.4 



Comparison 



0.2 0.4 0.6 0.8 1 

Frequency (mHz) 

Figure 3. Fourier transform of the J1630 light curve (top panel). 
The data shows low-frequency sky noise at ~0.2 mHz. After pre- 
whitening, the strongest peak is at the orbital period of the bi- 
nary (dashed line). No significant peak is observed at the same 
frequency for the comparison star. 
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Figure 4. High speed photometry of J 1630 over 9.6 hours (top 
panel). The bottom panel shows the same light curve binned into 
100 phase bins. The solid line shows a 0.17% amplitude signal at 
the orbital period. 



cal analysis, the doppler boosting signal is likely detected in 
J1630. However, the boosting signal does not provide any 
new physical constraints on the properties of this binary. 



4 DISCUSSION 

Along with J0106 and J0651, J1630 is only the third de- 
tached WD binary known to have a period shorter than an 
hour. All three were discovered in the last year as part of 
the ELM Survey. Previously, all known systems with P < 1 
hr were interacting AM CVn systems. We do not see any 
evidence of mass transfer (no emission lines and no obvious 
excess continuum) in the three targets mentioned above. The 
primary WDs in J0106 and J0651 both show ellipsoidal vari- 
ations due to their larger size and shorter orbital periods. 
These variations are extremely useful for constraining the 
inclination angle of the systems. However, J1630 does not 
show any significant ellipsoidal variations, and we only have 
an upper limit on the inclination angle due to the lack of 
eclipses. 

For the companion to avoid detection in the SDSS pho- 
tometry and our spectroscopy implies that it is ^ 10 x 
fainter than the visible WD. For an edge on orbit, such a 
companion woul d have M = 0.30 M ^, T e s ^ 7500 K, and 
R ^ 0.021 R l|Panei et alj 120071 ). A total eclipse would 
be ^70% deep and last for about 1.8 minutes. The lack of 
eclipses in the photometry constrain the inclination angle 
to i ^ 82°. Hence, J 1630 is best explained by a binary sys- 
tem containing a 0.30 M WD with a M ^0.30 M WD 
companion at a separation of ^0.32 Rq. 

The two WDs in the J1630 binary will merge 
in ^31 Myr due to gravitational wave radiation 
l|Landau fc Lifshitz I ll95SD . After J0651, J1630 is cur- 
rently the second quickest WD merger system known. 
When the mass transfer starts, J1630 will most likely 
have unstable mass transfer due to the mass ratio o f 
its components being close to unity ( M arsh et al.1 12004| ) . 
However, the merger outcome is uncertain because of the 
unknown inclination angle and the companion mass. If the 
companion is another He-core WD, the merger will likely 
create a single He-burning subdwarf in 23-31 Myr. If the 
companion is a more massive carbon/oxygen core WD, the 
system will merg e in gC23 Myr to form a rapidly rotating 
massive WD (see lKilic et al.| [2010, for other possibilites) . 

Short period binary WDs are important gravitational 
wave sources. For example, the 12 minute orbital period bi- 
nary J0651 sh ould be detected by LISA within its fir s t wee k 
of operation l|Brown et al.l l201ld ; iNelemans et al.1 120041 ) . 
Similarly, the previously discovered 39 min orbital period 
system J0106 may be detected by LISA after 1 yr of observa- 
tions. With an orbital period similar to J0106, J1630 is also a 
promising candidate for detection. For an average inclination 
angle of 60° and model-dependent distance of 0.7 kpc, we ex- 
pect the gravitational wave strain at Earth log h = — 22.0 at 
a frequency logv (Hz) = -3.08 l|Roelofs et al.l 120071 ) . This 
is at the S/N = 5 detection limit of LISA after 1 year of 
observations. 
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5 CONCLUSIONS 

We report the discovery of a new 40 min orbital period de- 
tached binary WD system, J1630. Along with the previously 
discovered 12 min and 39 min orbital period systems, J0651 
and J0106, J1630 is only the third known detached binary 
WD system with a period less than an hour. All three sys- 
tems are excellent gravitational wave sources to be detected 
by LISA. 

J1630 is almost a twin of the 39 min orbital period 
system J0106. However, due to its smaller size and larger 
mass, significant ellipsoidal variations are neither expected 
nor detected in its light curve. Therefore, the inclination 
angle of the system cannot be constrained accurately with 
the current data. Follow-up high-speed photometric obser- 
vations at a larger telescope will be useful to constrain the 
doppler boosting signal more accurately and to search for 
grazing eclipses. The two WDs in this system will merge in 
less than 31 Myr and likely form a single subdwarf star or a 
rapidly rotating massive WD. 
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